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Starting from an exact relationship between news, threshold and price return distributions in 
the stationary state, I discuss the ability of the Ghoulmie-Cont-Nadal model of traders to produce 
fat-tailed price returns. Under normal conditions, this model is not able to transform Gaussian news 
into fat-tailed price returns. When the variance of the news so small that only the players with zero 
threshold can possibly react to news, this model produces Levy-distributed price returns with a —1 
exponent. In the special case of super- linear price impact functions, fat-tailed returns are obtained 
from well-behaved news. 

PACS numbers: 

Few agent-based models of financial markets are amenable to mathematical analysis. Notable exceptions include 
ercolation-based imitations models , behavioural switching models 0, , and Minority Game-based models 0, 
Mathematical solutions are valuable tools for exploring the mechanisms underlying financial market dynamics, 
particularly in a view to understand the causes of unpredictability and large fluctuations. Thresholds mechanisms, 
whereby the agents react to events exceeding their personal tolerance levels, provide an indirect way to avoid forcing 
the agents to trade at each time step, and also implicitly define timescales. 

A new simple Markovian model of agents that adjust their thresholds on past volatilities was recently introduced 
and exhibited promising features 0. Here, I derive exact relationships between price return, thresholds and news 
distribution in the limit of infinite number of agents, and discuss in details when this model produces fat-tailed returns 
from well-behaved news distributions. 

The model works as follows |6j: there are N agents with individual threshold &i, i = 1, • • • , N, who monitor a single 
source of external news. At time t the news intensity and sign is e(t), drawn from the symmetric distribution p(e); a 
positive e corresponds to good news, and reversely. Only the agents whose thresholds are smaller than the absolute 
value of the news react to it by placing a order to buy if e(t) > and reversely. Mathematically, if <fii(t) denotes the 
action of agent i at time t, 

= e(e(t))0(e t - 0((t)) - e(-e(t))e(0 4 - e(f)), (1) 

where Q(x) is the Heaviside function; note that <fii(t) G { — 1,0,1}. This gives rise to global demand imbalance 
$(f) = J2 i (f>i(t). The price evolves according to 

logp(t+l)=logp(t) + r(t), (2) 

where r is called the price return. Its relationship with demand imbalance is assumed to be fully characterised by a 
constant price impact function g 

r{t) = g($(t)/N). (3) 

•i-H , 

At each time step, agent i may re-adjust his threshold 6i with probability s to \r(t)\. In that way, the agents are 
given minimal learning capabilities, tracking changes of volatility with a reaction time scale of order 1/s. It is worth 
noting that the agents have no influence on the news process. The inventors of this model produced evidence that 
the agents transform Gaussian news into price return distributions P(r)whose shape and properties depend on the 



parameters N , s and p, some of them seemingly not Gaussian anymore 
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Stationary state 

Since p{e) is symmetric, the model can be simplified by assuming that replacing e by e' = |e| and and p(e') by 
p'(e') = 2p(|e|). Dropping the primes, the dynamics of individual actions reduces to 



<Pi(t) = e(e t -0 i (i)). 



(4) 
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Because of Eq. (|3J), since $ is discrete for finite TV, r(t) and take discrete values. The master equation for P(8,t) 
is therefore 

P(0, t + 1) = (1 - s)P(6, t) + sS[6 - g( P(0\ *))] (5) 

8><u 

Letting t = 9/N and e' = e/N, and dropping the prime, one can transform this problem into a finite continuous 
interval (r, r and e belong to [0, 1]), resulting in 

P(t, t + 1) = (1 - s)P(t, t) + sS[t - g( dr'P(r', t))]. (6) 

Jo 

By definition, P(t, t) becomes P{t) in the stationary state; the final step is to average over all possible e, which yields 

PW - /**Wr-,iFMi) = P[F -fr-M)]'yb--(x)) - (7) 

It should be noted that the disappearance of s in this equation is a sign that s is a time scale; the auto-correlation 
function of r{t) stays positive for a time scale of order 1/s, as can be seen in the original paper. In the stationary 
state, since r = g(F(e)), the pdf of r has exactly the same functional form, because the thresholds are adjusted to the 
last price return with probability s. The multiple appearances of P(t) under various forms in the right hand side of 
Eq (J7J) is due to the self-referential nature of the model. 

The dynamics produces P(t) from p(e). However, Eq (0) seems quite difficult to solve; it is much easier to find 
which p(e) corresponds to a given P(r), since 

p(e) = P[g(F(e))]g'(F(e))P(e) (8) 

does not depend on itself in an intricate way. 



Linear price impact function 

The original model reports simulations with g(r) ~ Xr, which simplifies the discussion. Without loss of generality, 
one takes A = 1. The original model also assumes Gaussian news arrival, that is, p{e) = 2A/"(0, L>), e > 0. While it 
would hard to find P(r) that leads exactly to such a distribution, P(r) = 2A/"(0, D) corresponds to 

2 erf [vfe] 2 _ , 2 
L> z ir 

In particular, for e > D, p(e) approximates very well a Gaussian. It is therefore possible to conclude that Gaussian 
news produces Gaussian-tailed price return distributions. Numerical simulations show that Gaussian news can produce 
one peak or two symmetric peaks with or without one central Levy distributed peak (in finite size systems) whose 
origin is explained below. 

Nevertheless, one might question the assumption that news are Gaussian. For instance damages caused by natural 
disasters have fat tails 0- The power-law case yields a more direct relationship between news and price returns: 
P(r) = (a + X)r a (a > —1) corresponds to p(e) = (a + l) 2 e a ( Q + 2 ) = (a + l)e' 3 for r and e 6 [0, 1]: it is therefore 
impossible to produce a < from fj > with linear price impact functions, hence, well-behaved news cannot trigger 
fat-tailed price returns; however, fat-tailed news do trigger fat-tailed prices. Fig ^ confirms the validity of the above 
calculations: numerical and theoretical values of P(r) match perfectly. The case of P(r) is interesting: even after 
10 5 time steps, there is a discrepancy between numerical simulations and theoretical predictions; this comes from 
the temporal correlations of the thresholds over time scales of order 1/s, which causes the presence of peaks in any 
snapshot of P(r). When one averages P(r) over several snapshots, one finds that the theoretical prediction and the 
numerical simulations match better. Alternatively, one could plot the distribution of the values taken by the threshold 
of a given agent. As a final example, an exponential distribution of thresholds P(r) = e~ T , r > 0, corresponds to 
p(e) = e e e~ e . 

So far, the discussion has only addressed the model in the limit of infinite systems. However, the behaviour of the 
model changes completely for finite size if D <C 1/N and sN > 1. In that case, it is easy to convince oneself that the 
only situation yielding non-zero price returns occurs when P(0) ^ 0; in other words, in the case of very small news, the 
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FIG. 1: Threshold (left) and price return (right) cumulative distribution function F(t), respectively F(r) obtained with 
p(e) = (a + lfe a{a+2) for a = 1 (black symbol and line) and a = 2 (red symbol and line). N = 10000, s = 0.001, 100000 
iterations; F(r) is averaged over 9 equally spaced snapshots in the last half of the iterations. 




FIG. 2: Cumulative distribution of r and r for TV = 5000, D = 0.1/N, s = 0.1, 10 s iterations. 

only agents who react to them are infinitely susceptible, as they react in fact to any news amplitude. The mechanism 
is the following: the smallest non-zero threshold is 1/N > e* for all t. Assume that P(9 — 0,to) = 0, which implies 
that r(to) = 0. Then on average sN agents adjust their threshold to 8 = 0. At the next time step, all of them react 
to any news amplitude, which leads to an average price return r(to + 1) = s. Because of the threshold adjustment 
process, the number of agents with a zero threshold follows a multiplicative random walk r(t + 1) = r(t)rj(t) where 
V(rj) — 7V(1 — s, s(l — s)/N), with a barrier at 1/N; when P(0,t') < 1/N, r jumps back on average to s. Numerical 
simulations fully confirm this mechanism. They also show that P(r) a r _1 (and of course P(t) oc t -1 ), with an 
upper cut-off at about s. Multiplicative random walks with lower boundary are well known to give rise to power-law 
distributed density functions whose exponent 7 is given by (see for instance 0) 

J dr/P(? ? )?7 7 " 1 = 1. (10) 

In the large N limit, V{rj) converges to a Dirac function, which leads to (1 — s) 7_1 =1, hence 7=1. The normalisation 
of P(r) is not problematic, since r is bounded from above. Intermediate cases occur when sN ~ 1: the centre of the 
distribution is of Levy type with —1 exponent, and is truncated. 
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Non-linear price impact functions 

Average real market price impact functions are reportedly concave functions, either g(x) — x 1 with 7 < 1 |gL 1 1 0l | or 
g(x) = log a; 0). It will be enough to discuss the power-law case. Using P(r) — (a + l)r a (a > —1), one finds that 

p(e) = 7 (q + i)2 e 7(o+D 2 ~i = 7 ( a + !)2 e ^_ (n) 

Equivalently, a — — 1- The condition for power-law distributed price returns (a < 0) is (3 < 7 — 1. Therefore, 

sublinear impact functions dampen /?, and fat-tailed returns are only obtained from sufficiently fat-tailed news; indeed, 
if 7 — 1 < (3 < 0, the news are fat-tailed, but the price returns are not. A superlinear g on the other hand makes it 
possible to have a well-behaved p and fat-tailed price returns, but, as stated above, average price impact functions 
are not convex functions. 



Conclusion 



In short, in the N — > 00 limit, the above derivation not only shows the crucial importance of the functional form 
of the news on the price properties, but also that this model produces very fat-tailed price returns in the limit of 
infinitcsimally small news with exponent 1, and only for finite size systems; the other possibility for fat-tailed price 
returns arises when the price impact function is super-linear. Therefore, under realistic assumptions, this model 
cannot transform Gaussian news into fat-tailed price returns. 

I wish to thank Andreas Ipsen for stimulating discussions and an anonymous referee for useful suggestions. 
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